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Q Abstract 
< 

We discuss the proposal of Freedman, Ringo and Dombeck [jl] to search for 
the neutron electric dipole moment by use of the acceleration of ultracold 
neutrons in an inhomogeneous electric field followed by amplification of the 
, resulting displacement by several methods involving spin independent inter- 

actions (gravity) or reflection from curved (spin independent) mirrors. We 
show that the proposed technique is inferior to the usual methods based on 



o 



■ magnetic resonance. 

PACS: 21.20.Ky, 14.20.-c, 03.75.Dg 

V : I. INTRODUCTION 

73 . 

Since the appearance of our original paper [Q] there has been some additional discussion 
H and UJ], so that it now seems reasonable to publish a revised version with a more detailed 
exposition of our quantum mechanical treatment (sec. II B.) 

Searches for a neutron (or elementary particle) electric dipole moment (edm) are inter- 
esting because an observation of an edm would be a demonstration of the violation of time 
reversal (T) invariance outside the K° system. Freedman et al. [Ij] have proposed a new 
method to search for the neutron edm. The method is claimed to offer the possibility of 
vastly improved sensitivity due to the amplification of the effects of the interaction of an 
edm with an electric field by means of subsequent motion in a gravitational field or reflection 
from a convex mirror. In the present paper we will review this proposal and show that the 
claimed gain in sensitivity is based on a misunderstanding of a semi-classical model of the 
processes involved. 

According to the proposal, Ultra-cold neutrons (UCN) polarized along the x axis enter 
a chamber (accelerator) where the electric field directed along z has a gradient in the x 
direction. Since the incident spin state which is an eigenstate of a x can be considered as 
a coherent superposition of the two eigenstates \±) z of cr z , we can suppose that neutrons 
possessing a non-zero edm in one of the eigenstates will suffer an acceleration 
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After time T\ the two spin states will be separated by a distance 

After the period of acceleration, the neutrons are allowed to leave the chamber in a vertical 
direction (+z) rising against gravity. After reflection from a surface inclined at 45°, the 
horizontal separation Sx is converted into a vertical separation 5x = 8z. The neutrons then 
follow parabolic trajectories under the influence of gravity. Because of the difference in 
initial heights, 8z, hence a difference in kinetic energy, the two spins will accumulate a phase 
difference along the two trajectories 

A<f> = \ [ AVdt = mgAzT 2 (3) 



h 

where T 2 is the time of flight along the parabolic trajector(ies). This is called a 'gravitational 
amplifier' by the authors ]I]. It is then proposed to measure this amplified phase difference 
between the states |±) z as a precession of the polarization vector in the x, y plane: 

A0 = 9T 2 ^Tt (4) 

Putting in practical values for the parameters the authors expect a sensitivity to the edm 
of 10~ 28 e — cm, which is superior to that expected for other methods ||. It is seen that the 
proposed effect depends crucially on the description of the spin 1/2 system where the phase 
difference is calculated along trajectories ending at different points while the polarization is 
calculated by assuming the two states combine coherently at a single point. 

The authors also propose a second type of amplifier based on repeated reflections from 
a curved surface, given by Z(x). Then the two 'points' representing the two spin states 
separated by 5x due to the edm acceleration (|1]) will reflect from portions of the surface 
with slightly different slopes, differing by: 

d 2 Z (x) 

a W = dx 2 6X ( 5 ) 

and, after reflection the angle between the trajectories will increase by 2a (x). Again, the 
increased separation is supposed to result in an amplification of the phase difference and 
hence of the precession in the x, y plane. This has the same feature as the "gravitational 
amplifier": trajectories ending at increasingly distant points are used to calculate a phase 
difference which is then thought to be measured as a precession of the polarization which is 
calculated by considering that the two states combine coherently at a single point. 



II. DISCUSSION OF THE PROPOSAL ACCORDING TO DIFFERENT MODELS 

As the present proposal is based on a mixing of models we will present a discussion of 
the proposal from several view points. It is important that the models be well defined and 
applied in a consistent manner. Mixing different models leads to errors and confusion. 

A detailed discussion of the relation between classical, semi-classical and quantum me- 
chanical descriptions of similar situations has been given in || with a more complete quan- 
tum mechanical discussion in j7|. We will review these ideas with emphasis on the application 
to the present model. 
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A. Classical Model 



As an edm interacting with an electric field behaves identically to a magnetic moment 
interacting with a magnetic field we choose to discuss the problem in terms of magnetic 
moments in a constant magnetic field as this is probably more familiar. We will see later 
that the introduction of field gradients in the author's proposal is neither essential nor 
desirable. 

A classical magnetic moment jt coupled to an angular moment j '• , ~ft = 7 j in an 
external field B, (of course the same discussion will apply to an edm in an electric field) 
obeys the equation of motion 

d7 — -1? (6) 



dt 



fl x 



whose solution is seen to be a precession of j around B with the Larmor frequency u>l = jB 
independent of the angle between ~~f and B and which is a constant of the motion. If the 
field exists in a region of length L, then the particle will cross the field in a time t = L/v 
during which time the spin will precess through an angle 

ip L = v L t = u L L/v = ^BL/v (7) 

The spin components, a X)V , averaged over the velocity spectrum / (v) of the beam, will 

then be given by 

K> = J dvf (v) cos (8) 

V 



(a y ) = J dvf (v) sin ( 



This is the basis of a technique that is often used to measure velocity distributions. 

Note that as oj^L increases enough the polarization is expected to decrease. In a 
classical model the separation of states in space does not appear. As the gravitational inter- 
action is spin-independent there is no gravitational amplification in this model. Moments 
entering the field directed perpendicular to it undergo no energy change since jt ■ ~S = 0. 

The classical model is expected to give accurate results so long as the separation between 
trajectories associated with different states is small compared to the correlation (coherence) 
lengths of the wave function. Thus the classical model breaks down in the Stern-Gerlach 
effect but gives a very good description of Larmor precession. [0. Therefore one could 
conclude on this basis that the proposed amplification does not exist but we are aware that 
this argument would not be seen as compelling. 



B. Quantum mechanical treatment 

Since the semi-classical model is in some ways the most difficult, primarily because it is 
prone to misinterpretation, we will first discuss the quantum mechanical model. Note that 
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to our knowledge no quantum mechanical demonstration of the alleged amplification effect 
has been presented, ||, |4j]. 

We start with some general, elementary remarks. A wave function ip (x) and its Fourier 
transform \l/ (fc) are respectively, the probability amplitudes for finding the particle in a given 
region of position or momentum space. Then a displacement of the particle is equivalent to 
a phase shift (linear in momentum) of the momentum wave function: 



d 3 ke< 



k ■ af 



ip (l? + 5x 
A more general phase shift 

V (IT) = Itf fie) 



fic) 



k ■ 



(fx) 



up 



can be considered as leading to a displacement, 5x, such that 

ip (k) = ip (^k ) + 1? ■ ^k<P ( & ) = V 9 fi o) + ' & x 



(9) 



(10) 



(11) 



only under the condition that ^ ( k ) is narrow enough that higher order terms in the 

expansion of <p ( k ) can be neglected. Also if ( k ) should not play an important role. In 
the case of neutron spin echo || (or an edm accelerator with constant electric field) where 
(p oc l/k the condition for this is (n/k <C 1) where k is the width of the wave in momentum 
space centered around k Q . In the spin echo case we also have 5x ^> 1/ac so that the concept 
of a displacement between the states has physical meaning and, as will be seen below, no 
polarization can be observed under these conditions. It is necessary to cancel the phase 
shift (this is called obtaining the echo), and thus eliminate the displacement 5, in order to 
observe the polarization. 

The reason for belaboring these rather self-evident points is that in the work in question 
there is a tendency to make a distinction between the displacement and the precession meth- 
ods of searching for an edm whereas we have seen that the two concepts are only different 
ways of looking at the same phenomenon. There is only one phase for each quantum state 
and attempts to separate this phase into various components can often lead to confusion. 
This applies to the present case where a term linear in k is singled out as a 'displacement' 
or in attempts at separating the phase into 'geometric' and 'dynamic' parts. This latter 
separation is strongly dependent on the coordinate system used for the calculation - only 
the total phase remains unchanged ||. 

If we consider a wave function ( k j for a state where the spins are initially polarized 



in the x direction 



|*(t = 0)> 




(12) 



and, as a result of a spin dependent interaction, the spin states pick up an additional phase 
(p± so that at some time t we have 



$ (Te,t)) = * (Tc) 




(13) 
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then the expectation value of !? will be 

(a x +ia y ) = (*(t)\<r x + i<T v \*(t)) = J \v(t)\ 2 e<rt-^)d 3 k (14) 

In the case where the spin dependent interaction, V± = ±V , is small compared to the 
kinetic energy of the beam particles (this is certainly the case for any interaction involving 
a particle edm and holds for the spin echo case as well 0]) and slowly varying we can use 
the WKB approximation to write 

<p ± = f dx^Jk 2 - 2mV±/h 2 ^kxT I dx mV °Jf^ (15) 

J J fx K 

For V (x) = const = fiB or fi e E we can write, taking a = ULLm/2h, with lol = 2fiB/h or 
2n e E/h 

ft = ±r = ± (<Po ~ k6x) (16) 

where we have expanded around k Q , (k = k Q + k) for a narrow spectrum centered on k Q 
using Sx = a/kl, see (ITT]) , {ip Q = a/k Q ). The case where V a (x) has a constant gradient, 
V Q = n e x (ff ) is seen not to introduce any significant differences into fllSD . Then 



(a x + ia y ) = (x - 5x, t)tp(x + 8x, t) e i2ip °) (17) 
= (r+(x,t)^(x,t)) = I(t) (18) 

where ip (x) is the Fourier transform of \& (k) and equ. (|I7|) is the result of the Wiener- 
Khintchin theorem applied to (|14"D. We see that (a x + ia y ) in equ. (|14"D as a function of 
ip represents the Fourier transform of the beam velocity spectrum in agreement with the 
classical calculation, equ. (|j). From equ. ( |T7| ) it follows that this is a decreasing function of 
ifk = ft ~ fk or $ x f° r large enough ipk so that increasing <pk or Sx far enough will result in 
a reduction of the net polarization. As 5x increases so much that the correlation function in 
(|j~7l) approaches zero we approach the case of the Stern-Gerlach effect where the spin states 
can be considered as truly separated and one can talk about a displacement without any 
confusion. Thus we have seen that the discussion can be carried out equally either in terms 
of (fk or 8x, they both represent the same physical situation. 

The main thrust of the proposal under discussion |l[ is that after a period in the electric 
field region the x velocity of the UCN is changed into the z direction by reflection from a 
mirror at 45°. Subsequent motion in the gravitational field is supposed to significantly in- 
crease the phase difference (p, according to the (incorrect) semi-classical argument presented 
in M and outlined in the introduction. 

The key point is that no spin-independent interaction can influence the phase shift be- 
tween the two spin states. At the entrance to the gravitational field drift region (called the 
"amplifier" by the authors) the wave function will be of the form given by fllTf ) and (|i~6|); 
taking the Fourier transform 
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for a narrow spectrum. The superscript ° refers to the incoming beam at the entrance orifice 
of the amplifier region. 

We now calculate the wave function at the exit of this region. It is sufficient to consider 
only the motion in the z direction and confine ourselves to the steady state situation. Then 
using the WKB method for the case of the slowly varying gravitational potential, which is 
not necessarily small compared to the particle kinetic energy, we can write the wave function 
at the output (in the case of zero electric field in the accelerator) 

$(z) = ^(jfc^/o VC* 5 ^* ^A^e^WVA* 1 (20) 

k k 

where K 2 = 2m 2 gz/h 2 = m 2 v 2 /h 2 . v z is the velocity of a particle that falls a distance z, 
starting from rest. In order to investigate the relation between phase shifts and displacements 
we wish to study the behavior of ^ in a small region located around z = Z, so we set z = Z+e, 
e<Z and expand the exponential in (E^) around Z. 

e (21) 

Since A (k) represents a fairly narrow wave packet centered around k Q we can write 
fc = !;„ + k, k C fc < Kz (the latter condition is necessary for the 'amplification' to be 
significant, see below) and then we obtain from keeping only the terms in k and e 

ip (z = Z + e) « A (k + k) exp i 

K 

The terms in fl22|) can be interpreted as follows: The first term represents a displacement in 
position by v Q tz, the distance the particle with the initial velocity, v Q , would go in the time 
tz that the particle takes to fall to Z and an additional displacement E /mg associated with 
the boundary condition at z = 0. The next two terms represent the change in wavelength 
due to the acceleration during the fall. Since J2 K A (ho + K ) e%Ke represents the envelope of 
the wave function at z — 0, the last term in (^2|) shows that the envelope is spread by the 
factor 7] = Kz / k Q ^> 1 due to the energy dependence of the index of refraction (dispersion) 
for the Schrodinger wave. 

When the electric field is switched on we have (using equ. fllEp.) 

i)± (z = Z + e) w A ( k o + K ) e^ Kt e iip i 

K 
K 

displaying only the relevant terms. Now if we calculate the center of the wave packet, e , 
according to d<p/dn\ eo = 0, we find the center of the packet is indeed shifted by 

5e = r]5x (24) 

in agreement with the ideas of ref. The wave packet spreading by the factor rj, effects 
the electric field induced displacement by the same factor during the transit through the 
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ip{z 



A (k) exp i 



{k 2 + Klf 2 + (k 2 + K^ 



3m 2 g 



1/2 



h 2 k k 2 k 

^(K z -k )K + Kze + 7 ^e + ^Ke 
m A g 2Kz Kz 



(22) 



(23) 



amplifier region. This is the quantum equivalent of the semi-classical argument given in [|TJ . 
However we see immediately from ( Effi) that the k5x term in the phase shift is a small part 
of the total edm induced, spin dependent phase shift and neither term in the phase shift is 
altered by travel through the 'amplifier' region. Thus from equ. (^) we have 

_ « E A * ( k o + K ') A ( k o + «) e^ (K_K,) ' e W<Po-(«-«')Sx) (25) 

{z—Z+e) 

k,k' 

We see that the edm induced phase shift at the output is exactly the same as at the 
input to the 'amplifier'. Thus any measurement performed subsequently, whether using a 
polarization analyzer or interferometer with magnetic beam splitters or any other scheme, 
will yield the phase shift as it was at the input and there is no amplification. Any attempts 
to measure the increased displacement directly as a displacement are easily seen to be 
completely impractical as is evidently recognized by the authors of UJ. 

In the next section we show that a correct semi-classical argument leads to the same 
conclusion. 

In response to the circulation of an earlier version of this paper, |2| , Peshkin || introduces 
what he calls the "No-Go" theorem. Since if)± (jt,t) satisfy the same Schroedinger equation 
in the amplifier region (the Hamiltonian is spin independent) the quantity, I (t), (see equ. 

(0)) 

I (t) = J d 3 x t) tp- ( x>, t)] (26) 

is independent of time due to unitarity. 

However this leaves open the possibility that the integrand over some limited region 
of space might be time dependent, thus allowing the amplifier to work without violating 
the theorem. Thus Peshkin states that measuring the polarization "in a small range of 
x ..instead of over all x at one time , to avoid the integral over all space in I (t) ..avoids 
the No-Go theorem in principle" . He also claims that the No-Go theorem does not address 
interferometer experiments as "there the phase shift shows up as an overlap integral between 
two partial wave packets in one emergent beam only, not as the conserved integral over all 
space. For the same reason, the theorem also does not speak usefully to versions of the 
proposal |I| in which the phase shift is measured with a Mach-Zender interferometer instead 
of a polarimeter" . The same argument is presented in a recent Letter to the Editor of Nuc. 
Instr. and Meth. in Physics Research, by Dombeck and Ringo, PJ. The idea is that at 
the output of an interferometer (with the spin states separated by a magnetic mirror and 
travelling through the different arms of the interferometer) the output (considering that one 
spin state is flipped inside the interferometer) would be given by 



the cross terms giving the integral I(t), equ. ( pB"D but with the integral taken only over one 
of the emergent beams. However the partial beams ip± (~af , t) at the output of the interfer- 
ometer share all the same phase properties of the beams at the input of the interferometer 
and the theorem will apply equally to both situations. Be that as it may we propose a more 
stringent version of the 'No-Go' theorem, the 'Never-Go' theorem, i.e. the motion of the 
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spin is unaffected by a gravitational field due to the spin independence of the Hamiltonian 
as we have shown above (equs. ^ and |25]) 

We have shown that while an amplification of the 'displacement' betwen the two spin 
states does occur in the 'gravitational amplifer' due to the wavelength dependence of the 
index of refraction, this is virtually unobservable; the phase shift between the two spin states 
remains unchanged on traversing the region as does the direction of the neutron polarization. 

As we have not integrated over space or time our result shows that attempts to avoid 
the No-Go theorem by confining the measurements to limited regions will not work. The 
inapplicability of the theorem is not a sufficient condition for the existence of an amplification 
effect. 



C. Semi-classical model 



This model employs the geometrical optics approach to quantum mechanics, where each 
spin eigenstate is represented by a different trajectory. In a sense this is the most difficult 
model as it is prone to misunderstanding. In [1], at least three errors are made in appli- 
cation of the semi-classical model to the spin amplification process. Let us first consider 
amplification by vertical displacements, as discussed before equ. (3) above. 

First, assuming that the concept of a displacement of the two spin eigenstates is correct, 
we can calculate the gravitational effect on the spin precession angle. The first error in 
[1] occurs when the assumption that the phase between the two wave functions is simply 
the phase difference between the two eigenstates evaluated at the respective maxima of the 
wave function envelopes. This is a quantum mechanically incorrect procedure because the 
phase determination does not commute with the determination of the wave function center 
at a given time; in other words, it makes no more sense to compare the phases between the 
two eigenf unctions at two distinct spatial points than it does to compare the phases at two 
different times. 

The correct procedure is to make a point by point comparison between the two wave 
functions, then average over the two envelopes. This is the procedure normally used when 
calculating the interference between two scalar or vector fields as is commonly done in 
electrodynamics (see, for example, §|, Sec. 7.2). That this is the correct procedure can also 
be seen from the fact that detection occurs at a single point in space-time (for example, 
the neutron is absorbed on a 3 He nucleus thereby "collapsing" the wave function to a single 
space-time point; the spin direction is given by the two wave function phases at that point 
in space-time). 

We can now properly calculate the phase difference between the two eigenstates at a 
fixed point in the final polarimeter/detector; it is the phase difference between two classical 
trajectories that meet at the same space-time point, initially separated a distance 5x per- 
pendicular to the momentum k. The change in phase is simply the change in action along 
the two trajectories, and this can be easily calculated to first order by use of a theorem, 
which is of crucial importance to interferometry (but universally ignored) due to Chiu and 
Stodolsky flO |. This theorem states that a change in the action when one of the endpoints 



of a classical trajectory is displaced is given by 

SS = P?5x° - P^xl (27) 



S 



with summation notation over spatial coordinates implied (/i = x,y,z), and where S is 
the action (equal to the quantum mechanical phase up to a factor of h), refers to the 
momentum, x^ the path endpoint coordinates, refers to the trajectory beginning, and D 
the trajectory end (at the polarization analyzer /detector). As discussed already, Sx^, that 
is, the relative displacement of the path at the detector, is identically zero because a neutron 
is detected at a single point (e.g., a polarized 3 He nucleus). The displacement of the path 
starting point is given by 5x x = 5x as given by equ. (2) above. However, it is assumed in [1] 
that 5x is perpendicular to the neutron momentum; therefore, the change in action is exactly 
zero, as given by the Chiu and Stodolsky theorem, and there is no gravitational acceleration 
effect, in essence, by definition. Any other effects that could change S, particularly those 
relating to the electric field gradient or gravitational acceleration in the specific geometry 
given in [1], enter only in second or higher order. 

The above arguments can be immediately applied to the curved mirror amplifier. We 
again assume that the trajectory endpoints must meet in order for there to be an interference, 
and we calculated the change in action as above. Again we find that SS is identically zero 
to first order in the electric field gradient. 

The final semiclassical misconception in [1] concerns the use of an electric field gradient 
over the storage volume. It seems to us that a larger (at least two times) 5x (hence phase 
shift) can be generated by sending the "bipolarized" neutrons from a region of zero electric 
field to a region of constant high electric field. Each eigenfunction will acquire a change 
in energy, hence a change in velocity, as it enters the electric field region; although 5x 
only increases linearly in time, for a given storage time Ti, suddenly accelerating the two 
eigenstates to their final velocities would lead to a larger 5x than if the two eigenstates were 
subjected to a weaker electric field gradient averaged over the storage volume, but giving 
the same final velocities only after storing for a time 7\. The implication is that the use of 
an electric field gradient is completely pointless and only leads to a dilution of a possible 
edm effect. However, this is not surprising based on the foregoing considerations: to achieve 
the maximum sensitivity to an edm effect, whether that be interpreted as 5x or a change 
in phase between the two spin eigenfunctions, the interaction energy given by the usual 
Hamiltonian 

H = - t iJ'E (28) 

must be as large as possible over the duration of an experiment, for it is the integral of H over 
time that gives the relative action between the two eigenstates. We thus see immediately 
that if the average magnitude of E is compromised by the wasting of electric field strength 
toward the establishment of gradients in the system, the final net sensitivity to the edm 
interaction given by H must also be compromised. 



III. CONCLUSION 

We have analyzed the proposal for a new type of neutron edm experiment from three 
different perspectives, and in each case, have arrived at the same conclusion: The new tech- 
nique offers no gain in sensitivity as compared to the usual magnetic resonance technique. In 



9 



fact, a careful analysis reveals that the new method is inferior to the conventional methods. 
In [1], a semiclassical approach was improperly used to analyze the proposed technique, 
and a comparison with the atomic interferometer of Kasevich and Chu |TTJ was used to 
justify the approach. However, there really is no point of comparison between the atomic 
interferometer and the system described in [1]. The Kasevich and Chu "interferometer" is 
based on a superposition of internal quantum states of the Cs atom, specifically, the ground 
state hyperfine levels; there is no discussion here of the center of mass wave function, but 
only of the phase difference between the internal states. This phase difference evolves at 
the hyperfine frequency (approximately 10 GHz) and the beauty of the system is that the 
freely falling atom experiences a Doppler shift relative to an oscillator fixed in the labora- 
tory; this relative frequency shift in the accelerating system makes possible, for example, 
a precise measurement of the Earth's gravitational field. In a certain sense, the Kasevich 
and Chu system really isn't an interferometer (this point is addressed in [12|]); the evolving 
internal quantum state phase difference serves as a clock which can be compared to the 
stationary laboratory oscillator, and the system is best described by the classical approach 
given above. We might invoke the quantum or semiclassical model to calculate the result of 
some force that would cause the two hyperfine levels to spatially separate; the result of this 
calculation would simply show a diminishing of the internal interference effect because the 
two hyperfine eigenfunctions no longer fully overlap in space-time; in the limit where the 
separation is greater than the center of mass wave function coherence length, the concept of 
a superposition of internal quantum states entirely loses its meaning. 
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